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We present a theory of the Josephson effect in a twofold-tilted Josephson junction made by d-wave 
anisotropic layered superconductors. We find the appearance of an intrinsic electrical resistance that 
arises from the misalignment of the superconductive planes (the CuO^ planes in YBCO) in the two 
electrodes. This intrinsic contribution to the tunnel barrier has several non-trivial consequences. 
The result is relevant for understanding the electric transport properties of [100] tilt and [100] tilt-tilt 
Josephson junctions based on d-wave superconductors. 

PACS numbers: 74.50.+r, 74.72.-h, 74.78.-w 



Grain boundary Josephson junctions (GBJs) have been 
fabricated and studied immediately after the discovery 
of High Tc superconductors Nonetheless the mech- 
anism of the supercurrent transport in these devices is 
still matter of scientific investigation and argument of 
general physical interest [2], |3j. It is well known that 
the dc Josephson effect can be described by means of a 
scattering theory using the Bogoliubov-de Gennes (BdG) 
equations Q. Characteristic localized states of quasi- 
particles, known as Andreev bound states, are found in 
the superconductive energy gap region through which 
the Josephson supercurrent flows between the two elec- 
trodes [5]. This kind of theoretical approach has been 
successfully used in describing analytically the physics 
of in-plane grain boundary cuprate Josephson junctions 
([001] tilt) where a two-dimensional scattering theory 0- 
|8j finely incorporates both the propagation of the quasi- 
particles along the Cu-0 planes and the anisotropic d- 
wave-induced symmetry of the pair potential [§J . More- 
over in ref 2] a numerical approach, also based on the 
BdG equations, has enriched the understanding of the 
role of charge inhomogeneities in limiting supercurrents 
in [001] GBJs. 

However, a number of exp eriments have been carried 

out on [100] YBCO GBJs[lM3> where the Cu -° P lanes 
of one or both electrodes are tilted by an arbitrary angle 
((f>i and (f>2 in Fig. [IJ with respect the substrate surface, 
implying that the relevant geometry of the junction is 
no longer in-plane. In particular, relatively high values 
of the I c R n products [12Hl4| have been observed in such 
YBCO GBJs. [100] junctions seem to be very promis- 
ing for developing new sensors of radiation in the tera- 
hertz band 1161 Il7j | and observing macroscopic quantum 
properties [l5j . Indeed Kawabata et ai [l8| have pointed 
out that macroscopic quantum effects in YBCO devices 
may be better observed in junctions showing electrodes 
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with pair potential lobes aligned each other (do/do junc- 
tions). This configuration can be obtained with out-of- 
plane [100] GBJs geometries, whereas [001] GBJs are not 
suitable for fabricating do/^o junctions [9(. 

In spite of all these circumstances, as of today, and 
differently from the [001] tilt junction specific- 
theoretical treatment of the transport properties in [100] 
junctions is lacking 19]. One can probably individuate 
emblematically the obstacle of developing a model the- 
ory for this case in the non-trivial modification suffered 
by the quasiclassical trajectories of the excitations. In- 
deed such trajectories are no longer confined to a single 
plane as in [001] layered junctions 0, but lay in 

planes tilted with respect to the substrate plane by the 
misorientation angles <pi,4>2i see FigJT] 

In this Letter we take an important step towards the 
effective evaluation of the peculiar transport properties 
of YBCO [100] tilt junction (a L =a R = in Fig. [J) and 
[100] tilt-tilt junctions (junctions in which one or both 
angles oil, an are different from zero). 

We focus on the dc Josephson effect and find that the 
influence of the different tilting of the conduction planes 
with respect to the barrier plane (on passing from one 
electrode to the others) brings about an enhanced normal 
state electric resistance of the junction. The key result is 
an analytical expression for the Andreev spectrum for the 
quasiparticles, which retains the influence of the tilting 
of the Cu-0 planes as well as the effects of the d-wave 
anisotropic symmetry of the pair potential. 

In order to describe charge transport in a meso- 
scopic [100] YBCO grain boundary Josephson junction, 
we assume a superconductor-normal-insulator-normal- 
superconductor (SNINS) model @, HH-Sl]. The super- 
conductive material considered is a d-wave anisotropic 
layered cuprate. The presence of normal regions of the 
order ~ £o = hvp/irAo, the ballistic coherence length 
[24| . is introduced for modeling the mechanism of the 
Andreev reflection and a possible suppression of the or- 
der parameter near the junction. 

Furthermore we suppose that the quasiparticles are 



2 



<P R =<P 




FIG. 1: a) z axis view of a layered SNINS [100] tilt-tilt junc- 
tion. The thickness of the normal region is a+b, the d-wave 
superconducting regions occupy the space defined by x < —a 
and x > b. The insulating barrier (I) is located at x = 0. b) 
3d view of the conducting planes with the indicated tilt an- 
gles 4>\ and (j>2- The d-symmetry pair potential is represented 
by A L (6) = A cos(26» - 2a L ); A R (6) = e iv A cos(26> - 2a R ), 
f = fR — fL is the global phase difference between the two 
superconducting regions. The solid and dashed lines repre- 
sent the electron-like and the hole-like elementary excitation 
trajectories extending over the two planes, respectively. 



constrained to move exclusively along the Cu-0 planes 
of the two electrodes. This last assumption is based on 
the fact that the normal conductivity in c axis direc- 
tion in YBCO is about one hundred times smaller than 
that along the ab plane. In the superconducting state 
this strong anisotropy persists in the supercurrent distri- 
bution owing to a large ratio Xl/X 2 b where A c (A a &) is 
the London penetration depth across the planes (in the 
planes) (3. |23|. 

The junction barrier, the plane x = in the frame xyz 
indicated in Fig[TJ is normal to the substrate and is con- 
sidered perfectly flat. This reflects the fact that in a 'val- 
ley' type morphology junction, like the one represented 
in Fig. [T]d, there is almost perfect matching of the con- 
ducting planes at the GBJ [26]. These latter are tilted 
with respect the substrate plane z — 0, by the angles </>i 
and —(f> 2 in the left (x < 0) and in the right (x > 0) elec- 
trode, respectively. The two SN interfaces are normal to 
the x axis, as indicated in Fig. [TJ The coupled motion 
of hole-like (v) and electron-like (u) components of the 
wave function >F in the Cu-0 planes is described by the 



2d quasiclassical Bogoliubov de Gennes equations |27j ]. 
For the plane x' — y and x > they write (analogous 
equations hold in the plane x" — y for x < 0) 

Eu(x', y) = h(x', y)u(x', y) + A R (k)B(x' - b)v(x' , y) 

Ev(x',y) = -h(x>, y)v(x', y) + A R (k)Q(x' - b)u(x', y) 

(1) 

where h(x',y) = -h 2 (d 2 /dx' 2 +d 2 /dy 2 )/2m+V(x)-E F . 
0(2) is the Heaviside step function and Ep = H 2 kp/2m 
is the Fermi energy. The insulating barrier is mod- 
eled by a delta function, V(x) — Uq5(x), Uq be- 
ing the Hartree potential. The spatial dependence of 
the pair potential is assumed to have a step-functional 
form described by A(fc, r) = 0(— a — x)Ai(fc) + <d(x — 
b)An(k). The bound states solutions of Eqs (pj for 
and *(x") [^(x',y) = <F(x') exp (ik y y), *(x", y) = 
<F(x") exp (ik y y)] in the Sl, Nl, N r and Sr regions, in- 
cluding normal and Andreev reflections of holes and elec- 
trons, are respectively 
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P denoting right (R) or left (L), fc£' e , k± h are the wave 
numbers of the electron-like and hole-like quasiparticles 
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respectively, see FigQJi, that move in the superconduct- 
ing regions, kf, the wave numbers of the electrons and 
holes in the normal regions respectively, and E < | | . 
The conducting layers in the electrodes are treated as 
planes parallel to the x" — y plane or to the x' — y 
plane. The propagation coordinates x', x" in the left and 
right electrode are given by x' = x cos 0i + z sin 0i , x" = 
x cos 4>2 — z sin </>2 , respectively. For d-wave symmetry the 
effective pair potential is modeled as = Ao| cos(2(6*=f 

a L ))\e^± , A^ = A | cos(2(6»Tafl))|e i(v+ ^± ) [IQ, where 
ip± are the phases as felt in a normal reflection by the 
electrons (+) and holes (-), ip is the global phase differ- 
ence between the two superconducting regions and ccl, 
an are the angles between the crystallographic a-axis of 
the left and right superconductors and the normal to the 
y axis. 9 is the incident angle of the quasiparticle tra- 
jectories to the y axis. Boundary conditions have to be 
specified at the three different interfaces, x — —a, x = 
and x = b, to determine the twelve unknown coefficients 
in ^f(x') and These conditions are: the conti- 

nuity of the function 'J and its derivatives, <M> /dx' and 



d^/dx", across the clean interfaces x — —a, x — b, re- 
spectively. Moreover, at the normal-insulator-normal in- 
terface (plane x = 0) where the change in the propagation 
direction from x" to x' occurs, boundary conditions are 
the continuity of '5 and the discontinuity of d^! /dx im- 
posed by the presence of the (5-function barrier insulator. 
Then, since d/dx = cos 4>\d/dx' and d/dx = cos 02G? / dx" , 
we write the matching conditions at x = as 
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Using Eq. @ and imposing the above conditions, un- 
der the assumption of perfect retro-reflectivity {Ep ^> 
E, |A£|) of the Andreev reflections [2(J, we derive an ho- 
mogeneous linear system of equations for the the twelve 
unknown coefficients. The condition of existence of solu- 
tions, in addition to the trivial one, provides the following 
spectral equation: 
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and where the effective barrier parameter Z e ff, defining 
the electrical transparency of the junction T e ft = 1/(1 + 
Zl }f )), appears 



Z eff {9, 01, 02) 



Uz{9) 2 + (cos 0i -cos0 2 ) 

4 COS 01 COS 02 



(7) 



with Z(6) = Z /cos9 and Z a = kpUo/2Ep the usual 
barrier parameter. With the definition Eq. ([7]) we may 
write the spectral equation in the same form which holds 
for the [001] case described in [2(| and (23|. However Eq. 
([?]) shows a new contribution to the barrier, coming from 
the specific geometrical configuration depending on the 
angles 0i and 02- Equation (|5|) provides the quasipar- 
ticle energy levels as a function of the superconducting 
phase difference <p for the case of layered SNINS d-wave 
[100] tilt-tilt junctions and other limits (SIS, SNS, INS). 
When the junction is in the clean limit, (Zq = 0) the 
transparency of [100] junctions does not reach the uni- 
tary value but it is still limited to 

4 cos 01 COS 02 



1 eff 



(COS 01 



which is 6- independent. Eq. ([8]) embodies the pure effect 
on the junction transparency of the misalignment of the 
conducting planes between the two electrodes. The zero- 
temperature normal-state conductance G n is evidently 
also affected by the angles 0i and 02 and can be written, 
in terms of the average transparency (T e ft), as 



^{Teff) 



(Teff) 



1 eff 



C 2 coth- 1 (VI + C 2 ) 



(9) 



where C — 2Zo/(cos0i + cos 02), and 2(T e //) = 
/-ir/2 d9 cos{9)T e f f(9 : 0i, 02). Therefore, the effect of the 
tilting of the planes in SNINS junction manifests as an 
excess of electrical resistance in the normal state. This 
effect persists also in the limit of zero size normal regions 
( a = b = 0). 

One of the consequences is that the Andreev level spec- 
tra are strongly modified. Indeed, in FigJ5] the Andreev 
spectra for two different GBJ configurations are reported. 
In particular we have considered, most representatively, 
the case of [100] tilt-tilt mirror (q^ = —cir = a) GBJ, 
with 0i = 0, 02 = 7r/4, and a = tt/8 (dashed line) 
and the analogous mirror [001] junction (solid line) both 
in the clean limit (Z = 0). As is well known, in mirror 
junctions, for a given symmetric rotation a (0 < a < 45°) 
around the c-axes, depending on the quasiparticle inci- 
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FIG. 2: Comparison between Andreev level spectra of d-wave 
SNS 45° mirror (ol = —or = 22.5°) grain boundary junc- 
tions calculated in the clean limit (Zq = 0) for [001] and [100] 
configurations, respectively. Solid and dashed curves refer 
to [001] GBJs with fa = 2 = 0, and [100] GBJs (0i = 
and 02 = 45°) respectively. Curves a): edgegap-like state 
Andreev levels for a trajectory with quasiparticle incidence 
angle 8 = 7r/10; curves b) midgap-like state Andreev levels 
with 6 = 3tt/16. 
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FIG. 3: Comparison of the temperature dependence of the 
absolute value of the Josephson critical current in d-wave [001] 
tilt (solid curves) and [100] tilt-tilt mirror (oll = — ur = a) 
grain boundary junctions with barrier parameter Zo = 0.7 
for three representative misorientation angles a. a) a = 13°, 
b) a = 18°, and c) a = 21°. For a — 18° arrow indicates 
the sign inversion of the critical current (0 — n crossover) in 
[001] GBJs. The remaining curves in b) show a shift of the 
crossover temperature with the tilt angle 02 in [100] GBJs. 
The sign of the critical current 7 C is positive for curves a), 
negative for curves c). 



dence angle 9(—tt/2 < 9 < tt/2), two kinds of bound 
levels exist [2211281]: midgap-like states, for ±7r/4— |a| < 
6 < ±7r/4 + |a| and edgegap-like states in the comple- 
mentary intervals. The formation of midgap states, i.e. 
the formation of zero energy states when the phase differ- 
ence across the junction is zero, is characteristic of d-wave 
Josephson junction 29]. Both kinds of levels determine 
the magnitude of the Josephson current as well as the de- 
pendence of the maximum Josephson current on the tem- 
perature 0- The two kinds of Andreev levels, selected 
for two illustrative 6 values, have been derived according 
to Eq. [S] In particular, when the levels are degenerate 
in energy, i.e. at <p = n, an asymmetric tilt of [100] junc- 
tions (0i 7^ 02) opens up an energy gap (E gap ) in both 
kinds of Andreev level spectra similar to what happens 
in the presence of impurities. For selected trajectories 
(6 = and 9 = ir/4), the entity of this bandgap may be 

^gap/An = 2t/l 



simply expressed as E„ 



T e ° // |sin2a| 



for 



i; ^ | cos 2a | 



edge-like energy states and Eg ap /A = 2 

for midgap-like states, respectively. 

The Josephson current may be derived directly from 
the energy spectrum. As a consequence it is sensitive 
to the modifications derived so far. In the short junc- 
tion limit, for which a + b << £ Uj the discrete Andreev 
levels determine all the Josephson current through the 
expression [si |30| 
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where the index n labels the bound Andreev energy lev- 
els and f(E n ,9) is the Fermi function. As a matter of 
fact it turns out that the Josephson current shows a de- 
pendence on the angles 0i and 02 in [100] geometries. 
Fig. [3] shows a comparison of the dependence on the 
temperature of the maximum Josephson current between 
asymmetric [100] and [001] mirror junctions. The cur- 
rent is normalized with respect to Iq = 2ekpL y lS.Q/h, 
which is the zero temperature maximum Josephson cur- 
rent through an s-wave SNS junction in the clean limit 
[30l | . The superconductive gaps in the two electrodes 
are assumed to obey a Bardeen-Cooper-Schrieffer (BCS) 
temperature dependence. It is well known that in this 
kind of junction the interplay between midgap and edge- 
gap states may lead to a temperature sign inversion of 
the maximum Josephson current (Fig. [3] b)), i.e. to a 
temperature dependent — 7r crossover [7J, |22| (arrow in 
Fig. [3]). For the chosen value of the barrier strength pa- 
rameter Zo = 0.7, the modifications for increasing values 
of the tilt angle 02 are evident, showing a shift towards 
higher temperatures of the — 7r crossover as 02 increases 
from 0° to 45°, see curves b) in Fig. [3] 

In this Letter we have presented a theoretical model 
for the description of electrical transport in YBCO [100] 
tilt valley type GBJs based on the quasiclassical BdG 
equations. We have compared our results with the well 
studied [001] tilt case. We have derived the spectral rela- 
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tion in the case of [100] tilt-tilt junctions. The Andreev 
bound states calculated through this equation show mod- 
ifications analogous to those caused by the presence of 
an insulating layer. However the nature of this modifi- 
cations is different: it depends on the geometrical mis- 
match between conducting planes. We have evaluated 
the temperature dependence of the maximum Josephson 
current in the case of [100] mirror junctions and found 
a dependence of the 0-to-7r crossover temperature on the 



inter-electrode misorientation angles. 
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